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1 . SCOPE 


This document contains the design specifications and description 
of the modifications made to the Statistics and Data Transformation 
processors of the EOD-LARSYS system. The requirements to be sat- 
isfied are specified in the task agreement 77-3 titled, "LARSYS/ 
Data Transformation Processor Modification'* originally dated 6/17/77 
and revised 9/12/77 and the "EOD-LARSYS System Design Modifica- 
tions for PCG Transformation Requirements" memo, dated 8/24/77. 

The modification requirements were provided by the Research, Test, 
and Evaluation Branch (RT&E) of the Earth Observations Division 
(EOD) , of NASA/ JSC. 



2 . APPLICABLE DOCUMENTS 


The following documents will serve as references; 

• EOD-LARSYS Users' Document, LEC-3984, March 1977 

• Task Description and Agreement, dated 6/17/77, "LARSYS/Data 
Transformation Processor Modification" 

• Symat, Covar, Test Procedures for Matrix Calculations by 

W. L. Morris, C. L. Wiginton, and D. K. Lowell, Report #17, 
Contract NAS 9-12777, dated October 1972 

• Principal Component Greenness (PCG) Transformation Requirements 
for LARSYS by R. A. Abotteen 

• EOD-LARSYS System Design Modifications for PCG Transformation 
Requirements, dated August 24, 1977 


3 . SYSTEM DESCRIPTION 


This document describes the changes that will be implemented 
into the existing Statistics and Data Transformation processors 
to generate and output an optional greenness and/or Principal 
Component Greenness (PCG) Image. The basic image will continue 
to be a multi-temporal/multi-pass data tape. The result is a 
transformation which will assist in LACIE-type classification. 

The user will be provided the options of : 

a. Filtering raw or sun angle corrected data vectors in the 
Statistics processor, 

b. Using sun angle corrected data in the DATA - Transformation 
processor to compute a greenness and/or principal component 
image . 

3 .1 HARDWARE DESCRIPTION 
N/A 

3.2 SOFTWARE DESCRIPTION 

3.2.1 SOFTWARE COMPONENT NO. 1 (STAT) 

The Statistics processor, STAT, will be modified to calculate on 
option the mean and covariance matrices from filtered and/or sun 
angle corrected data. These options will be provided if the new 
control cards 'SUNANG' and 'OPTION WATER' or 'OPTION TASSEL' are 
input by the user. 

If the 'SUNANG' control card is present, the channel-oriented sun 
angle bias and gain correction factors will be applied to the cor- 
responding channel component of the raw data vector. 

Xy(I) = SAGAIN(I) * Xj^(I) + SABIAS(I) 

where X^(I) = the I-th channel component of the sun angle cor- 
rected data 

SAGAIN(I) = the sun angle gain correction factor for the 
I-th channel 



SABIAS(I) = the sun angle bias correction for the I-th 
channel 

Xj^(I) ** the I-th channel component of the raw data 

The user will have the option of applying one of the following 
filters to the raw or sun argl*^^ corrected data vectors: 

a. Water filter - initiated by the input of the 'OPTION WATER' 
control card. (see Appendix A. 2 for the 'Water* filter test) 

b. Tassel Cap filter - initiated by the input of the 'OPTION 
TASSEL' control card. (see Appendix A. 2 for the 'Tassel 
Cap' filter test) 

3. 2. 1.1 Linkages 

The STAT processor uses the FORTRAN-V compiler, Univac software 
system routines, EOD-LARSYS utility routines, and the C0MBK4 and 
STBASE common blocks. 

See Appendix A for the modifications to be made to the individ- 
ual subroutines, 

3. 2. 1.2 Interfaces 

The STAT processor interfaces with the EOD-LARSYS executive routine, 
MONTOR. 


3 . 2 . 1 . 3 Inputs 

The processor requires an MSS data tape (DATAPE) . The default 
tape assignment is the logical unit C (Fortran Unit 3) . 



• NEW CONTROL CARDS 


Keyword 
(col. 1) 
SUNANG 

SUNANG 


B~I*1ATR1X 


B-MATRIX 


THRESHOLD 


Parameter 
(col. 11-72) 

Tape 

(Default - no sun angle 
correction applied) 

n. r . . • ^ n . 

j 

(Default- no sun angle 
correction applied 

Cards 

(Default - none) 


Pile 

(Default - none) 


Sli/ • • • $ cl» 

1 1 
(Default - for 'Tassel 

Cap' filter, 

a^ “ 100 

a2 = 8.0 

a^ = 6.0 

a. - 10.0 
4 

a^ = 35.0 

for the 'Water' filter, 
a^ = 43.0 
- 12 . 0 ) 



Function 

Sun angles will be ex- 
tracted from the ERIPS 
'unload' MSS image tape. 

s sun angle for pass 
i, l<i<8. Sun angles are 
integer numbers which 
are 5<n^<85. 

The B-Matrix will be 
found on the cards im- 
mediately following this 
control card. 

The subroutine BMFIL 
will be called to read 
the B-Matrix off of the 
file, BMFILE 

The a^'s are decimal 
(floating point) numbers, 
separated by commas; if 
the user requested a 
'Tassel Cap' filter, j 
must be equal to 5; if 
the 'Water' filter was 
requested, j must be 
equal to 2. 


Keyword 

BIAS 


OPTION 

OPTION 


P arameter 

t • • . f 

(Default - b. » 0.0) 

X 


Tassel 

(Default - none) 


Function 

The b^'s are decimal 
(floating point) numbers, 
separated by commas; the 
bias values will be used 
only if the 'Tassel Cap' 
filter is requested. 

Initiates the Tassel Cap 
filtering. 


Water Initiates the Water 

(Default - none) filtering. 


^ MODIFIED CONTROL CARD 

Keyword Parameter 

(col. 1) (col. 11-72) 

CHANNELS Data = 1^^ , . . ., Ij^, 

(Default - n - 30) , 

Filter = m. , . . . m . 

I' . D 

(Default - none) 


Function 


The 1. 's and m. 's are 
1 1 

integer channel numbers, 
separated by commas; they 
must be a subset of the 
MSS Image tape. The Ij^'s 
will be used in the gen- 
eration of the mean and 
covariance matrix. The 
mj^'s are the channels to 
be used in the user-re- 
quested filter process. 

If the "Tassel Cap' filter 
was requested, j must be 
equal to 4 (based on 4 
channels/pass) . If the 
'Water' filter was requested 
j must be- equal to 2 . 


> 4 ^ 


3. 2. 1.4 Outputs 

The Statistics processor will continue to output the statistics 
on the default or user-requested SAVTAP file. The default logical 
unit is A (Fortran unit 1) and the default file number is 1, 

The line printer and file output remain unchanged. 

3. 2. 1.5 storage Requirements 
TDB 

3. 2. 1.6 Description 

The Statistics processor, STAT, of the EOD-LARSYS system will be 
modified to optionally calculate and output on SAVTAP a mean 
vector and mixture covariance matrix calculated from the op- 
tionally filtered raw or sun angle corrected data vectors. Several 
new control cards will be implemented in the STAT processor to 
determine if the above calculations are to take place. The pre- 
sence or absence of the 'SUNANG* control card will determine if 
raw or sun angle corrected data vectors are to be used in the genera- 
tion of the above matrices. The user will also have the option 
of filtering out some data vectors from the above calculations 
by the presence or absence of the 'OPTION TASSEL' or 'OPTION 
WATER' control card. 

The new mixture covariance matrix will be used in the DATA-TR 
processor to generate the mixture covariance matrix of the Green 
Image bands . 

3. 2.1.7 Flowchart 
N/A 

3. 2. 1.8 Listing 



TBD 


3.2.2 SOFTWARE COMPONENT NO. 2 (DATA-TE) 

The Data Transformation processor, DATA-TR, of the EOD-LARSVS 
system will be modified to perform an optional Green Image and/or 
Principal Component Greenness (PCG) Image transformation (s) . 

The addition of a SUNANG control card allows for the sun angle 
gain and bias correction factors to be applied to the raw data 
before the above mentioned transformation (s) . The Green image 
and/or PCG image will be additional to the current options for 
performing transformations in DATA-TR. 

3. 2. 2.1 Linkages 

The DATA-TR processor uses the Portran-V compiler, Univac soft- 
ware system routines, EOD-LARSYS utility routines, and the INFORM, 
TRBLCK, and GLOBAL common ?/iocks. 

3. 2. 2. 2 Interfaces 

The DATA-TR processor interfaces with the EOD-LARSYS executive 
routine, MONTOR. 

3. 2. 2. 3 Inputs 

The processor requires an MSS data tape (DATAPE) . The default 
tape assignment is the logical unit C (Fortran Unit 3) . 


• NEW CONTROL 

CARDS 



Keyword 
,(col. ,1) 

Parameter 
(col. 11-72) 

Function 


OPTION 

Green 

(Default - Green image 
not generated for out- 
put) 

The Green 
generated 

image will be 
and output. 



Keyword 
(col. 1) 

Parameter 
(col. 11-72) 

Function 

OPTION 

PCG 

(Default - No Princi- 
pal Component Green 
Image output) 

The PCG image will be 
generated and output. 

FORMAT 

Universal 

(Default - Universal) 

The transformed data 
will be output in the 
Universal format. 

FORMAT 

LARSyS 

(Default - Universal) 

The Transformed data 
will be output in the 
LARSYS II format. 

SUNANG 

Tape 

(Default - no sun 
angle correction 
applied) 

Sun angles will be ex- 
tracted from the ERIPS 
"unload" MSS image tape. 

SUNANG 

n^^/ . , t f Uj 

(Default - no sun 
angle correction 
applied) 

n. = sun angle for pass 
i, l<i<8. Sun angles are 
integer numbers » 5<nj^<85. 

TRPORM 

OUTPUT/UNIT*N, 

FILE=M 

(Default - N=14, 
M=l) 

N = the tape ’ s Fortran 
logical unit number that 
the new image of trans- 
formed data is to be out- 
put on. 

M “ the file number that 
the new image is to be 
output on. 



If both Greenness and PCG 
are to be generated, the 
former image will be on file 
M and the latter image on 
file (M+1) 


3. 2, 2. 4 Outputa 


The Data-Transformation processor will continue to output the 
transformed/rescaled data on the default or user-requested TRFORM 
file and/or unit, in either the Universal or hhRSYS II format. 

The default logical unit is L (Fortran unit 14) and the default 
file number is 1. 

The two additional optional transformer?/rescaled data output 
are Green Image and Principal component Greenness (PCG) Image. 

3. 2. 2. 5 Storage Requirements 
TBD 

3. 2. 2. 6 Description 

In order to allow the user to optionally generate and output 
a raw or sun angle corrected Green Image and/or PCG Image, 
several existing subroutines will be modified. In addition, 
several new subroutines and formulas will also be implemented 
into the Data Transformation processor. 

The scaling parameters MAX, MIN, and CON calculated in KBTRAN 
will be derived from the sun angle corrected SAVTAP subclasses ' 
mean and standard deviation terms when the sun angle correction 
option is requested. 

The new subroutine GETGCV will be called to pull out the mixture 

covariance matrix ^ from SAVTAP. The new subroutine PCSUNP 

X 

will provide the sun angle correction parameters. 

The green image will be generated as follows; 

G = Ax + b 

where G = transformed green image vector 

X = raw (or sun angle corrected) data vector 


A = matrix in which each row contains the LI, L2, or user- 
specified Kauth greenness vector (stored internally in 
BMAT in the 'B-Matrix' format) 

b = default or user-specified bias vector 

The covariance matrix will be used to generate the Green 

Image's covariance matrix according to the formula: 

where £„ = the Green Image's covariance matrix 

V5 

= the (optionally) filtered raw or sun angle cor- 
rected covariance matrix calculated in STAT 
T 

A = the transpose of the A Transformation matrix 

The newly generated covariance matrix S q will then be input into 
the new subroutine PCMMAT to derive an nxn unitary matrix PCM. 

The PCG image will then be generated as follows: 

P = (PCM) G 

where ? = n-dimensional PCG image vector (n - number of passes; 
i.e., LANDSAT acquisitions) 

PCM = an nxn unitary matrix obtained from PCMMAT 

G = n-dimensional green image vector 

The currently available options for re-scaling the transformed 
data to a range of 0-255 for output will be unchanged. All re- 
scaling options will be available to be applied to the green image 
and/or the PCG image created as a result of these modifications. 

3 . 2 . 2 . 7 Flowchart 
N/A 





3. 2. 2. 8 Listing 
TPD 
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A.l SOFTWARE SUBPROGRAM NO. 1 (SETUPl. ) 


The modifications to the subroutine SETUPl will result in the 
reading and decoding of the following additional control cards; 

a . SUNANG 

b. B-MATRIX 

c. THRESHOLD 

d. BIAS 

Two new additional parameters will be read and decoded from the 
'OPTION' control card. They are: 

a . OPTION TASSEL 

b. OPTION WATER 

The format of the 'CHANNELS' control card will be modified to be; 

CHANNELS DATA = d^, . . . , d^, FILTER = f^, . . . , f^ 

where the d ^ * s and f ' s are integers . 

Several new flags will be initialized in the subroutine to be 
tested in the LEARN subroutine. The flags that will be placed 
in the STBASE common block are; 

a. SAKEY - sun angle key - if on, apply the sun angle correction 

factors to the current raw data. 

b. FKEY - filter key - if FKEY=0, no filtering will be 

applied to the data vectors; if FKEY=1, the Water 
filter will be applied to the data vectors; if FKEY=2, 
the Tassel Cap filter will be applied to the data 
vectors . 

The following dimensioned variables will be initialized in this 
subroutine to be used in the LEARN subroutine ; 

a. SSUNAG - the sun angles read from the 'SUNANG' control card 

b. BIAS - the bias vector to be used in the Tassel Cap filtering 



process . 


c. TH - the threshold values to be used in the^Tassel Cap ^ 

or "Water’ filtering process. 

d. BMAT - the matrix (in the B-MATRIX format) to be used in 

the Tassel Cap filter test. 

e. FCIIN - the channels of the data vector that will be used 

in the user-requested filter test. 

A. 1.1 Linkages 

SETUPl is called by the Statistics processor driver program, 

3TAT . 

The following subroutines may be called by SETUPl; NXTCHR, 

FIND, NUMBER, BMFIL, and FLTMUM. 

A. 1.2 Interfaces 

SETUPl references the following common blocks; STBASE and GLOBAL. 
A. 1 . 3 Inputs 

The calling sequence to SETUPl remains unchanged. 


• NEW CONTROL 
Keyword 

CARDS 

Parameter 


(col. 1) 

(col. 11-72) 

Function 

SUNANG 

Tape 

(Default: no sun angle 
correction applied) 

Sun angles will be ex- 
tracted from the ERIPS 
unload MSS image tape . 

SUNANG 

n^^, . • . , Uj 

(Default: no sun angle 
correction app],ied) 

n^^ = sun angle for pass 
i, l<i<8. Sun angles are 
integer numbers which are 
5<n^<85 . 



Function 


Keyword 
(col. 1) 

B-MATRIX 


B-MATRIX 


THRESHOLD 


BIAS 


OPTION 


OPTION 


Parameter 
(col, 11-73) 

Cards 

(Default: none) 


File 

(Default: none) 


'’l' • • • » 

(Default: for ’Tassel 

Cap' filter, 

a^ = 100. 

b .2 ” 8.0 

a^ 6.0 

a^ = 10.0 

= 35.0 

For the 'Water' filter, 
a^^ = 43.0 
a2 = 12.0 

, . . . , 

(Default: = 0.0) 


Tassel 

(Default; none) 
Water 

(Default; none) 



The B-Matrix will be 
found on the cards im- 
mediately following 
this control card. 

The subroutine BMFIL 
will be called to read 
the B-Matrix off of the 
file, BMFILE. 

The a^'s are decimal 
(floating point) numbers, 
separated by commas; if 
the user requested a 
'Tassel Cap' filter, j 
must be equal to 5,* if 
the Water filter was re- 
quested, j must be equal 
to 2. 


The b^'s are decimal 
(floating point) numbers, 
separated by commas; the 
bias values will be used 
only if the 'Tassel Cap' 
filter is requested. 

Sets the flag (FKEY=2) to 
apply the 'Tassel Cap' 
filter to the data vectors. 

Set the flag (FKEY=1) to 
apply the 'Water' filter to 
the data vectors . 


• MODIFIED CONTROL CARD 


The following control card will replace the current 'CHANNELS' 
format . 


Keyword 
(col. I) 

CHANNELS 


MSS Image tape. The 1^'s 
will be used in the gen- 
eration of the mean and 
Covariance matrix. The 
mj^'s are the channels to 
be used in the user - re- 
quested filter process. 

If the 'Tassel Cap' fil- 
ter was requested, j must 
be equal to 4 (based on 
4 channels/pass) . If the 
'Water' filter was re- 
quested, j must be equal 
to 2 . 


Parameter 
(col. 11-72) 

Data = • . . » 

(Default: n » 30), 
Filter = m^, . . . 
(Default: none) 




m. 


Function 

The l.'s and m. *s are 

X Iv 

integer channel numbers, 
separated by commas; they 
must be a subset of the 


A. 1 . 4 Outputs 

Two new additional line printer outputs will be implemented into 
the subroutine SETUPl ; 

a. If the Tassel Cap filter was requested but the number of input 
filter channels is not equal to four or the number of input thres- 
hold values is not five, the following error message will be out- 
put: 

'The number of input filter channels or threshold values are 

I 

not compatible with the filter option requested-Terminate 
Execution ' . 


b. If the 'Water' filter was requested but the number of input 
filter channels or threshold values is not two, the above error 
message will be output. 

c. If the 'Tassel Gap' filter was requested but the user did not 
input the 'B-MATRIX' control card or if the dimension of the B- 
Matrix (KAUTH filter matrix) is not 4x4 then the following error 
message will be output; 

'The B-Matrix was not input or was not of the right dimension - 
Terminate . ' 

A. 1.5 Storage Requirements 
TBD 


A. 1.6 Description 

The purpose of the SETUPl subroutine is to read and analyze all 
of the input processor control cards as well as their parameters. 

As a result of the modifications to SETUPl, the user will be 
given the options to J 

a. Apply the sun angle correction factors to the raw data vectors^ 

b. Apply a 'Tassel Cap', 'Water', or no filter to the raw or 
sun angle corrected data vectors , and 

c. Calculate and output the filtered raw or sun angle corrected 
mean and covariance matrix computed over the image data set. 

The covariance matrix will be used in the DATA-TR processor to gen- 
erate a Principal Component Greenness Image. 

A. 1.7 Flowchart 
N/A 

A.l . 8 Listing 



TBD 


A, 2 SOFTWARE SUBPROGRAM NO. 2 (LEARN) 


The functions to be provided by the subroutine LEARN are to enable 
the user to calculate the mean vector and covariance matrix from 
filtered raw or sun angle corrected data vectors. The new flag 
SAKEY will be tested to determine if the sun angle correction 
factors are to be applied to the raw data vectors. The value of 
the new flag FKEY will determineif Water, Tassel Cap, or no filter 
is to be applied to the data vectors. 

A. 2.1 Linkages 

LEARN is called by STAT, the driver program for the STATISTICS 
processing. 

The following subroutines may be called by LEARN: TAPHDR, PCSUNF, 

PSBSFL, CMERR, LAREAD, LINERD, FDLINT, MATVEC, FLDCQV, FLDSPC, 
PLDHIS, CLSCOV, CLSSPC, CLSHIS , MULSPC, and SETMRG. 

A. 2. 2 Interfaces 

LEARN references the STBASE, GLOBAL, COMBK4 and ISOLNK common 
blocks. 

A. 2. 3 Inputs 

The calling sequence to LEARN remains unchanged. 

A. 2 . 4 Outputs 

There is no new line printer output from LEARN. 

A. 2. 5 Storage Requirements 



TBD 


A. 2.6 Description 

The modifications to be made to the subroutine LEARN are to test 
the sun angle (SAKEY) and filter <PKEY) flags. 

If the SAKEY flag is on (SAKEY » 1) , the following will take 
place: 

a. Test the flag ISUNT. If ISUNT is on (ISUNT « 1), the sun 
angles that relate to the B-Matrix channels will be read from 
the data tape. If ISUNT is off (ISUNT * O) , the sun angle 
were read in SETUPl. The sun angles will be used in the sun 
angle correction process. 

b. The subroutine PCSUNP will be called to output the sun angle 
bias and gain correction factors. 

c. The sun angle correction factors (bias and gain) will be 
applied to the raw data vectors before calculating the mean and 
covariance matrix. 

If the filter flag (FKEY) is equal to: 

a. Zero - the raw or sun angle corrected data vectors will not 
be filtered before calculating the mean and covariance matrices. 

b. One - the 'Water* filter will be applied to the raw or sun 
angle corrected data vectors before calculating the mean an6 
covariance. 

The data vector will be filtered if one of the following holds: 


where X. = the radiance value of the f. component of the current 

X X 

data vector, 
til 

= the i — channel value read from the FILTER parameter 
section of the CHANNELS control card. 



* the i~ user-supplied (the parameter read off of 
the THRESHOLD control card) or default (t^^ ■ 43.0, 
tg * 12.0) threshold value. 

c. Two - the Tassel Cap filter will be applied to the raw or 
sun angle corrected data vectors before calculating the mean 
and covariance matrices. 

The following transformation will be applied to each data vector 
before it is filtered: 



where Y = transformed raw or sun angle corrected data vector, 

K - the matrix input via cards or via the file, BMFILE, 

X » The components of the raw or sun angle corrected 

data vector, and 

b = the user - supplied or default (b^ * 0.0) additive 
bias vector. 

The data vector will be filtered if one of the following holds: 

1. tj^ 

2. Y„< t„ 

3. 

4. V 

5. Y4> tg 

The default threshold values are = 100.0, t2 = 
t^ = 10.0, and tg = 35.0. 


8.0, tg = 6.0, 


A. 2,1 Flowchart 
N/A 


A. 2, 8 Listing 
TBD 
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B.l SOFTWARE SUBPROGRAM NO. 1 (DATATR) 


The driver program for the Data Transformation processing, DATATR, 
will be modified to optionally generate an nxn unitary matrix PCM 
when the PCGC flag is on. The PCM matrix will be used in the 
generation of the Principal Component Greenness (PCG) Image, 

B.1.1 Linkages 

DATATR is called by the LARSYS executive program, MONTOR. 

The following subroutines may be called by DATATR: SETUPS, 

SETREM, KBTRAN, MAXMAT, TRHIST, LNTRAN, MTMLS6, MTMDAT, and 
the new subroutines GETGCV and PCMMAT. 

B.l. 2 Interfaces 

DATATR references the INFORM, TRBLCK, and GLOBAL common blocks. 

B . 1 . 3 Inputs 

The calling sequence to DATATR remains unchanged. 

B.l. 4 Output 

If the column dimension of the 'B-MATRIX' is not equal to the 
row and column dimensions of the data covariance matrix X, 

X 

read from SAVTAP, the following error message will be output: 

•THE DIMENSION OP BMAT XXX AND DCOV XXX ARE NOT COMPATIBLE - 
CALL EXIT ’ . 

B.l. 5 Storage Requirements 
TBD 


B.l. 6 Description 


The following modifications will take place only if the PCGC 
flag is on (i.e., the OPTION PCG control card was read in SETUPS) ; 



a. The subroutine GETGCV will pull the data covariance 

matrix SAVTAP. This covariance matrix will have 

X 

been calculated from filtered raw or sun angle corrected data 
vectors depending on the 'OPTION* and 'SUNANG' control cards input 
by the user during the STAT processor mode. 

b. Before the mixture covariance matrix of the Green Image 
bands ( 2 q) to be generated, the row and column dimensions 
(the row and the column dimensions should be equal) will be 
checked against the column dimension of the B-MATRIX, If these 
dimensions are not equivalent, the error message stated in the 
Output section 3. 2. 1.4 will be printed. 

c. The subroutines MTMLS6 and MTMDAT will generate the GCOV 
matrix. 

^G = B 2 

where ^ G ~ niixture covariance matrix of the Green Image 
bands, 

B = Transformation matrix in which each row contains 
the greenness vector. The matrix is to be input 
in the 'B-MATRIX' format, 

V - The optionally filtered data covariance matrix 
to be used in the PCG transformation will be 
read off of SAVTAP, and 

T 

B = the transpose of B. 

d. The subroutine PCMMAT will be called to operate on the matrix 
53o* The program will in turn return an nxn unitary matrix PCM. 

The PCG Image will then be generated by applying the PCM matrix 
to the Green Image. 


Q7 


B.1.7 Flowchart 
N/A 


B.1.8 Listing 


TBD 



B.2 SOFTWARE SUBPROGRAM NO. 2 (SETUPS) 


The modifications to subroutine SETUPS will result in the reading 
and decoding of the following additional and modified control 
cards: 

a. TRFORM OUTPUT/UNIT=n, FILE=m 

b. SUNANG n-j^f n2*..nj 

c . OPTION GREEN 

d . OPTION PCG 

The following flags will be set and used in the other subprograms. 
They will be added to the labeled common block, TRBLCK. 

a. PCF - flag indicating 'OPTION GREEN 'and/or 'OPTION PCG' 
control card was read in. 

b. SAKEY - sun angle key - if on, apply the sun angle gain, 

SAGAIN and sun angle bias, SABIAS correction factors ob- 
tained from the subroutine PCSUNF. 

c. PCGC - flag indicating 'OPTION PCG' control card was read in. 

d. GIC - flag indicating that the 'OPTION GREEN' control card 
was read in. 

The following variables will be added to the common blocks 
TRBLCK and GLOBAL to be referenced by the other subprograms; 

a. SAGAIN - sun angle gain correction factor. 

b. SABIAS - sun angle bias correction factor. 

c. TRNSFL - the file number that the output transformed data is 
to be output on (default = 1) . 

B.2.1 Linkages 

SETUPS is called by the driver program for Data Transformation, 
DATATR. SETUPS in turn calls the NXTCHR, FIND, BMFIL, NUMBER, 
ORDER, FLTNUM, CRDSTA, and WRTBM routines. If sun angle correc- 


tions for the Green Image and/or PCG Image Transformations were 
requested, a call will be made to the new subroutine, PCSUNF. 

PCSUNF will return the set of sun angle gain and bias correction 
factors to be used on each data channels input. 

B.2.2 Interfaces 

SETUPS references the following common blocks: INFORM, GLOBAL, 
ISOLNK, and TRBLCK. 


B.2.3 Inputs 

The calling sequence to SETUPS remains unchanged. 
• NEW CONTROL CARDS 


Keyword 
(col. 1) 

OPTION 


OPTION 


FORMAT 


FORMAT 


Parameter 
(col. 11-72) 

Green 

(Default - Green 
Image not generated 
for output) 

PCG 

(Default - no Principal 
Components Green Image 
output) 

UNIVERSAL 

(Default ~ Universal) 
LARSYS 

(Default - Universal) 


Function 

The Green image will be 
generated and output. 

The PCG Image will be 
generated and output. 


The transformed data will 
be output in the Universal 
format. 

The transformed data will 
be output in the LARSYS 
II format. 



Keyword 

Parameter 


(col. 1) 

(col. 11-72) 

Function 

SUNANG 

Tape 

(Default - no sun angle 
correction applied) 

Sun angles will be ex- 
tracted from the ERIPS 
"unload" MSS image tape. 

SUNANG 

, . . . , n j 

(Default - no sun angle 
correction applied 

nj^^sun angle for pass i, 
l<i<8. Sun angles are 
integer numbers, 5^n^<85 

TRFORM 

OUTPUT/UNIT=N, 

FILE=M 

(Default - N=14 
M=l) 

N = the tape’s Fortran 
logical unit number that 
the new image of trans- 
formed data is to be out 


put on. 

M = the file number that 
the new image is to be 
output on. 

If both Greenness and PCG 
are requested, the former 
image will be output on 
file M and the latter image 
on file (M+1) . 


B. 2,4 Outputs 

Modification to SETUPS will include the addition of the following 
error message; 

a. If the B-MATRIX input flag BMTRIG is not on (~1) , the error 
message will be: 

•The B-MATRIX CONTROL CARD WAS NOT FOUND - TERMINATE EXECUTION'. 

The line printer output will continue to list the (new) control 
cards read in by SETUPS . 


B.2.5 Storage Requirements 


TBD 


B.2.6 Description 

The function of SETUPS is to read and decode all of the input 
processor control cards as well as their parameters. The fol- 
lowing modifications will be made: 

a. The user will be able to specify on the new 'TRBORM' control 
card which unit number (default = 14) and file number (default = 
are to be used for the output transformed image tape. 

b. The user will be given the option of using raw or sun angle 
corrected data vectors in the transformation process depending 
on the presence or absence of the SUNANG control card. If the 
SUNANG control card was read in, the input 'MAXPT' vector will 
be sun angle corrected. 

c. The user has the option of creating a Green Image and/or 
Principal Component Greenness (PCG) Image. 

d. The transformed data will be output in the Universal format 
unless the LARSYS parameter is found on the FORMAT control card. 

e. If the Green Image and/or PCG Image is to be generated, the 
maximum dimension allowed for B-MATRIX is (n,n*k) , where k = 
number of channels per pass and n = the number of passes and 
n^*k<400. The default number of channels per pass (NCHPAS) is 
equal to 4. NCHPAS is set in the BLKCOM subroutine and included 
in the GLOBAL common block. 

B.2.7 Flowchart 
N/A 


B . 2 . 8 Listing 


B.3 SOFTWARE SUBPROGRAM NO. 3 (KBTRAH) 


The subroutine KBTRAN will be modified to test the new flag.' 

a. SAKEY - when it is on (SAKEY = 1), the channel-oriented 
sun angle correction bias and gain factors (obtained from the 
subroutine PCSUNF) will be applied to the respective channel 
component of the mean vector and covariance matrix obtained 
from the subclasses statistics file, SAVTAP. 

B.3.1 Linkages 

Subprogram KBTRAN is called by the Data Transformation driver 
program, DATATR wnen rescaling by the statistical method is 
requested by the user. KBTRAN calls the following subroutines; 
MATVEC, MTMLS6, MTMDAT, and PRTCOV. 

B.3.2 Interfaces 

The common blocks referenced by the KBTRAN subprogram are: 
INFORM, TRBLCK, GLOBAL, and C0MBK4 . 

B. 3 . 3 Inputs 

The calling sequence for KBTRAN remains unchanged. 

B. 3 . 4 Outputs 
If the SAKEY is on; 

a. The SAVTAP 's mean vectors will be transformed using the sun 
angle correction factors, SAG^IM and SABIAS. The formula used 
is : 

TM(I) 

where TM(I) 

SAGA IN (I) 

M(I) 

-'35 


= SAGAIN(I) * M(I) + SABIAS (I), 

= Ith component (channel) of the sun angle 
corrected transformed mean vector, 

= Ith sun angle gain factor that relates to the 
Ith channel, 

= Ith component (channel) of the mean vector 
from SAVTAP, and 


SABIAS(I) » Ith sun angle bias factor that relates to 
the Ith channel. 

b. The appropriate sun angle correction factors will be applied 

to the SAVTAP’s covariance matrix, , before transforming it 

V' T * 

to B , where B is the 'B-MATRIX' which contains the Kauth 

* T 

Greenness Vectors, and B is the transpose of B. The formula 
used on each component of £ is; 


TCOV(I,J) = SAGAIN(I) * SAGAIN(J) * COV(I,J) 

where TCOV(I,J) = sun angle corrected transformed (I,J)th 

component of the SAVTAP's covariance matrix 
I,J = 1, . . . , NOPEAT, 

SAGAIN(I) = Ith sun angle gain factor that relates to 
channel 1. 


SAGAIN(J) = J^th sun angle gain factor that relates to 
channel J, and 

COV(I,J) = (i,j)th component of the SAVTAP's covariance 
matrix. 


B.3.5 Storage Requirements 
TBD 


B.3.6 Description 

If the SAKEY flag is on, the sun angle correction bias and gain 
factors will be applied to the corresponding component of the 
SAVTAP mean vector and covariance matrix. 

If the PCGC flag is on (i.e. the 'OPTION PCG' control card was 
read in SETUPS) , the matrix PCM (output from the PCMMAT sub- 
routine) will be applied to the calculated Green Image's MAX^ 
(=MAX) , MINq(=MIN), and CONq(=CON) vectors to obtain the PCG 


Image's scaling parameters, MAXp, MINp, and CONp. 

The PCG Image's 

a. MAXp vector (PCM * MAX^) will be stored in locations BMAT 
(401) BMAT (401 + LCOMB-1) . 

b. MINp vector (PCM* MIN^) will be stored in locations BMAT 
(401 -f LCOMB) BMAT (401 + 2 * LCOMB-1). 

c. CONp vector (PCM * CON^) will be stored in locations BMAT 
(401 4- 2 * LCOMB) BMAT (401 + 3 * LCOMB-1) . 


B . 3 , 7 Flowchar t 
N/A 


B.3.8 Listing 


TBD 


B.4 SOFTWARE SUBPROGRAM NO. 4 (TRHIST) 


The following modifications to TRHIST will be performed if the 

PCGC flag is on (i.e. the 'OPTION PCG* control card was read in 

SETUPS) ; 

a. The scaling parameters, MAX, MIN, and CON, obtained from the 
histogram of the transformed data will be multiplied by the 
PCM matrix. The resulting vectors will be the PCG Image’s 
scaling parameters, AMAX, AMIN, and ACON. 

b. The PCG Image's scaling parameters will be stored in BMAT's 

scratch storage locations: BMAT (401) BMAT (401 + 3 ★ LCOMB-1) , 

B. 4 .1 Linkages 

The driver program DATATR calls TRHIST when the user specifies 

rescaling by the histogram method. 

TRHIST calls the following subroutines: LAREAD, TAPHDR, SQRT, 

PLDINT, LINERD, FDLINT, TRANSF, and MATVEC. 

B , 4 . 2 Interfaces 

TRHIST references the following common blocks; INFORM, TRBLCK, 

COMBK4, and GLOBAL. 

B. 4 . 3 Inputs 

The calling sequence to TRHIST remains unchanged. 


B.4. 4 Output ■ 

There is no line printer output. 
B.4. 5 Storage Requirements 



TBD 


B.4.6 Description 

TRHIST will obtain the scaling parameters, MAX, MIN, and CON for 
the transformed data from the histogram of the transformed data of 
the first field of a set of fields. 

If the PCGC flag is on (i.e. the 'OPTION PCG' control card was 
read in SETUPS), the matrix PCM (output from the subroutine 
PCMMAT) will be applied to the Green Image's scaling parameters 
MAXq (»MAX) , MIN (=MIN) , and CON^ (»CON) . 

The PCG Image's scaling parameters; 

a. MAXp (PCM * MAXq) will be stored in locations BMAT (401)”>- 
BMAT (401 + LCOMB-1) . 

b. MlNp (PCM * MINq) will be stored in locations BMAT (401 + 
LCOMB) '►BMAT (401 + 2 * LCOMB-1). 

c. CONp (PCM * CONq) will be stored in locations BMAT (401+ 

2 * LCOMB) BMAT (401 + 3 * LCOMB-1) . 

B.4.7 Flowchart 
N/A 

B. 4.8 Listing 
TBD 



B. 


B.5 SOFTWARE SUBRPOGRAM NO. 5 (PCMMAT) 


The new subroutine PCMMAT was originally progranuned for execution 
under EXECS by W. L. Morris, C. L. Wiginton, and D. K. Lowell 
(University of Houston Mathematics Department) . The revised 
subroutine will be incorporated in the Data Transformation sub- 
routine PCMMAT, PCMMAT operates on a real symmetric matrix A 
to produce an orthogonal matrix of approximate eigenvectors of A. 
For the Principal Component Greeness (PCG) transformation, the 
symmetric matrix will be the mixture covariance matrix of the 
green bands ( 2 q) and the output orthogonal matrix will be PCM. 
The n X n (n = LCOMB « row dimension of Z! q) unitary matrix PCM 
will have its ith row correspond to the ith largest normalized 
eigenvector of £q. The resulting components of PCM will be 
approximate eigenvalues (i.e. the error bounds will be between 
+ Q.OOOQQS) of In the present context, 

Eg = B Ey 

where E ~ mixture covariance matrix of the green 

bands . 

B = BMAT » matrix in which each row contains the 
greenness vector. 

y 

^ » the filtered raw or sun angle corrected 

data covariance matrix calculated in STAT, 

B.5.1 Linkages 

PCMMAT will be called by DATATR. 

PCMMAT will reference the following internal subroutines; MINDEX, 
ORDER, SUPSUM, and MATMUL. 

B.5. 2 Interfaces 

PCMMAT will not reference any common blocks. The interface with 
the calling program is by means of the calling sequence. 



B . 5 . 3 Inputs 

The callings sequence to PCMMAT is: 
Call PCMMAT (A, PCM, X, R, N) 
Parameter Dimension In/Out 

A (N,N) IN 


PCM (N,N) OUT 

C N OUT 

R N OUT 

N 1 IN 


Definition 

A symmetric real matrix 
(for PCG transformation 
processing, A=S the 
mixture covariance matrix 
of the green image bands) . 

An unitary orthogonal 
matrix of normalized 
eigenvectors of A. 

Vectors of centers of 
Weinstein discs. 

Vector of radii of 
Weinstein discs. 

The row and column 

dimension of A (for 
V 

La g , N=LCOMB=number of 
LANDS AT Passes) . 


B.5.4 Outputs 

The orthogonal unitary matrix PCM will be returned to the calling 
subroutine DATATR, via calling argument. 

P . 5 . 5 Storage Rcquircmonts 
TBD 


B . 5 . 6 Description 

proceeds through an iterative algorithm to produce the 



PCMMAT 


output orthogonal matrix PCM. PCM, which is composed of approxi- 
mate eigenvectors of q will call the following internal sub- 
routines : 

a. MINDEX - selects the order of operations within PCMMAT. 

b. ORDER - reorders' the components of the input vector into 
ascending order. 

c. SUPSUM - add the components of the reordered vector from 
ORDER. 

d. MATMUL - computes the matrix products. 

The relative error allowed in the approximate eigenvalues of 
Q will be set to 0.000005. 

B.5.7 Flowchart 
TBD 

B.5.8 Listing 


TBD 


B . 6 SOFTWARE SUBPROGRAM NO. b (TRANSF) 


The calling sequence to TRANSF will be modified to exclude the XT 
variable. TRANSF will be changed from a subroutine to a function. 
The internal variable XT will be set equal to TRANSF. The flag 
SAKEY will be tested to determine if the input raw data vector, 
IDATA is to be sun angle corrected before performing th ^ata 
transformation calculation. 

B* 6. 1 Linkages 

TRANSF is called by either the subroutines TRHIST and/or LNTRAN, 
TRANSF does not call any subroutine. 

B.6.2 Interfaces 


TRANSF refers to the TRBLCK common 
B.6.3 Inputs 

The calling sequence to TRANSF is; 

block. 

• XT(J)=TRANSF (BMAT, 

IDATA, TOP, 

IL, K, LCOMB, NS AMP, BIAS). 

Parameter 

Dimension 

In/Out 

Definition 

BMAT 

480 

IN 

The matrix, which is 
stored in the 'B-MATRIX' 
format, to be used in 
the transformation process. 

IDATA 

TOP 

IN 

The input raw data vector 
to be transformed. 

TOP 

1 

IN 

The maximum usable lo- 
cation in the array IDATA. 

IL 

1 

IN 

The component of the 
transformed data vector 
that will be generated. 



Parameter Dimension In/Out 

K 1 IN 

LCOMB 1 IN 


NSAMP 1 IN 

BIAS 16 IN 


Definition 

The current sample pixel 
being processed. 

The number of linear 
combinations or passes 
to be used in the calcu- 
lation. 

The number of sample data 
vectors to be transformed. 

The bias vector to be 
used in the data-transfor- 
mation process. 


B. 6. 4 Outputs 

There is no line printer output. The value of the IL-th compo- 
nent of the transformed data vector (XT)' will be returned in 
the function name TRANSF. 


B.6.5 Storage Requirements 
TBD 


B.6.6 Description 

The flag SAKEY will be tested in the function TRANSF to determine 
if the input raw data vector is to be sun angle corrected or not. 
If the flag was turned on (i.e. SAKEY = 1), the following trans- 
formation will take place before the existing data transformation 
process : 

TDATA = SAGAIN(I) * IDATA(L) + SAB IAS (I) 

TDATA = the sun angle corrected I-th channel com- 
ponent of the current data vector, 



where 


$AGA1N(I)« the sun angle gain correction factor for 
the I-th channel, 

IDATA{L)= the 1-th channel component ^ nsaSp+X ^ 
of the current data vector, 

NSAMP =s the number of sample data vectors 
(for the current scan line) that needs to 
be transformed, and 

SABIAS(1)« the sun angle bias correction factor for 
the I-th channel. 

If the SAKE! was not on (i.e. SAKEY = 0) , the raw data vector, 
IDATA will be used in the data transformation formula; 

XT ~ IDATA * BMAT + BIAS 

If the SAKEY was turned on (i.e. SAKEY = 1) , the sun angle cor- 
rected data vector, IDATA will be used in the transformation 
formula : 

XT = TDATA *BMAT + BIAS 

Before returning to the calling subroutine, TRANSF will be set 
equal to XT. 

B.6.7 Flowchart 
N/A 

B . 6 . 8 Li sti ng 
TBD 



B . 7 SOFTWARE SUBPROGRAM NO. 7 (LNTRAN) 


The modifications to be made to the sxibroutine LNTRAN will enable 
the user to generate and output a 'Green' and or 'Principal Com- 
ponent Greenness' (PCG) Image. The decision on whether these 
images are to be generated and output will depend on the new 
flags PCGC (turned on when the OPTION PCG control card was read 
in SETUPS) and GIC (turned on when the OPTION GREEN control card 
was read in SETUPS) . 

7. 1 Linkages 

LNTRAN may call the following subroutines: TAPHDR, LAREAD, 

PLDINT, WRTHDR, LINERD, FDLINT, TRANSP, WRTLIN, NTRAN, COMHST, 
and MATVEC. 

LNTRAN is called by the driver program of Data Transformation, 
DATATR. 


B.7.2 Interfaces 

LNTRAN will make reference to the following common blocks: INFORM, 

TRBLCK, COMBK4, and GLOBAL. 

B . 7 . 3 Inputs 

The calling sequence to LNTRAN remains unchanged. 

B. 7 . 4 Outputs 

The user will be able to generate and output a Green and/or PCG 
Image depending on the value of the flags GIC and PCGC. 

The transformed data set will be output on the default (Fortran 
Unit 14, file=l) or user-specified tape unit and file number. 

B.7.5 Storage Requirements 


TBD 


B.7.6 Description 


If the GIC flag is on (GIC=1 ) , the scaling parameters MAX, MIN, 
and CON derived from the histogram of the 'Green Image* will be 
used in the transforming, rescaling, and histogramming processes 
of the data vectors. 

If the PCGC flag is on (PCGC=1) , 

a. The scaling parameters (derived from the histogram of 'PCG' 
Image) stored in locations BMAT (401) BMAT (400 + 3 * LCOMB) 
will be used in the processes of transforming, rescaling, 

and histogramming the data vectors. Recall that the MAX 
vector is stored in locations BMAT (401) ->■ BMAT (400 + LCOMB) . 
The MIN vector in locations BMAT (401 + LCOMB) BMAT (400 + 2 
* LCOMB), and the CON vector in locations BMAT (401 + 2 * LCOMB) 
BMAT (401 + 3 * LCOMB) . 

b. The matrix PCM will be applied to the transformed data vectors 
output from the subroutine TRANST. 

B.7.7 Flowchart 
N/A 

B. 7 , 8 Listing 
TBD 



